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Abstract 

The simplest Higgs-portal dark matter model, in which a real scalar singlet is added 
to the standard model, has been comprehensively revisited, by taking into acconnt the 
constraints from pertnrbativity, electroweak vacnnm stability in the early Universe, dark 
matter direct detection, and Higgs invisible decay at the LHC. We show that the resonant 
mass region is totally exclnded and the high mass region is rednced to a narrow window 
1.1 TeV < ms < 2.55 TeV, which is slightly rednced to 1.1 TeV < mg < 2.0 TeV if 
the pertnrbativity is farther imposed. This high mass region can be fnlly detected by the 
Xenon IT experiment. 
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1 Introduction 


The Standard Model (SM) as the effective theory below the electroweak (EW) scale has 
been established after the discovery of SM Higgs scalar with mass around 125 GeV BE] 
at the Large Hadron Collider (LHC). Unfortunately, there is no viable candidate for a 
dark matter (DM) in the SM, which implies that an extension beyond the SM is necessary. 

Among viable extensions the simplest Higgs-portal dark matter (HDM) [3], IH [5l | 6 l [7] is 
the most economic and of special interest from both the DM- and LHC-phenomenology. 
In this extension, only a real singlet scalar s is added to the SM. With the aid of Z 2 
parity which is needed for the stability of singlet DM, the numner of model parameters is 
reduced to three. It implies that the minimal model is very predictable. In the subsequent 
studies on the HDM [H [9l [101 HB EM [IB EM [IB [IB [IB [IB [2B El] , the parameter 
space is analyzed in terms of experimental constraints as follows. 

• Dark matter relic density [22|, which shows that the low mass region, resonant mass 
region, and high mass region are all viable. 

• Dark matter direct detection such as XenonlOO [23], XenonlT [23] and LUX [25] 
experiments, which are able to effectively detect the high mass region up to ~ 3 
TeV. 

• Limit on the Higgs invisible decay h ss, which can effectively detect the low mass 
region with rus of order a few GeV. 

• Indirect detection such as limit on 7 -ray lines from the Fermi-LAT ^7] is rather 
efficient to constrain the resonant mass region [2B EB [30] . 

However, these constraints are not sufficient, as some possibly strong theoretic ones 
have been ignored. The object of this paper is to take into account all known experimental 
and theoretic constraints. Simultaneously, experiments constraints will be updated, if 
available. In compared with the previous works, e.g., a comprehensive review in [21], our 
main differences are the following ones. 

• The SM EW vacuum suffers from large quantum fluctuations in the early Universe 
as induced by the inflation, which imposes strong constraint [3B EB El] on the 
model parameters. SM vacuum stability must be taken into account to constrain 
the parameter space. For an earlier attempt, see [35l [36] - 
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• Assume the HDM as an effective theory between the EW and Plank scale, perturba¬ 
tive analysis should be considered seriously, which leads to strong upper bounds on 


the quartic coupling constant of I the DM self-coupling constant. This 


effective theory can be applied, e.g., to S-inffation, where the scalar singlet DM is 
identified as the inffaton m- 

• Very recently, the ATLAS collaboration has reported the latest limit about the 
Higgs invisible decay width |2B] in the low and resonant mass regions. This will be 
updated in our analysis. 

The paper is organized as follows. In Sec.2, we introduce notation and conventions in 
the HDM. In this section the constraint on the model parameters from the requirement 
of perturbativity is studied. In Sec.3, we discuss the constraint on the model parameters 
from the SM vacuum stability for 125 GeV Higgs mass. No additional mechanism for 
stabilizing the SM vacuum stability is employed except the DM scalar. 

In Sec.4 we re-examine the DM phenomenology. In particular, we calculate the s 
scalar relic density and the spin-independent nucleon-DM scattering cross section by using 
MicrOMEGAs [38], which are consistent with results in some earlier literature. In Sec.5, 
we update the LHG phenomenology about scalar DM, by focusing on the latest limit on 
the Higgs invisible decay. 

Finally, we present our main conclusions in Sec.6. We find that the resonant mass 
region is totally excluded and the high mass region is reduced to a narrow window 1.1 TeV 
< nis <2.4 TeV. If the requirement of perturbativity is further imposed, the allowed 
mass range is slightly reduced to 1.1 TeV < mg < 2.0 TeV. In either case the mass 
window can be fully detected by the Xenon IT experiment. 


2 HDM within Perturbative Region 


The Lagrangian for the simplest HDM with a Z 2 parity, under which s is odd and SM 
particles are even, is given by. 


^ — -^SM + 2 (^'®) + ^ 


( 2 . 1 ) 


where 



( 2 . 2 ) 
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In Eq. fl2.2p the first term denotes the Higgs potential with EW scale v ~ 246 GeV, and the 
A^-term and K^-ferm refers to s self-interaction and Higgs-DM interaction, respectively. 
Expand the helds as if = {v + h)/y/2 and s = (s)4-s = 0-|-s along the vacuum structure 
for positive and KsS, one obtains from Eq. fl2.2p . 


V{s,h) 


^ 2 2 I 1 I 

-m,s +ys + 




^-S^h+^S^h\ 
2 4 


(2.3) 


where rus = j2. Remarkably, there are only three new parameters in the HDM. 

It implies that this model can be studied concretely. 

In some situations peturbativity should be imposed on the model. This is true when 
it is assumed to the effective theory below some high energy scale such as Plank mass 
scale. See, e.g. [37] for an application of such idea to the S-inflation. 

In Fig{T]we show the upper bounds on the EW values of Kg and A^ by the requirement 
of perturbativity 


0 < < \/4^, 0 < As(/i) < \/4^, 


(2.4) 


for Renormalization Group (RG) scale /i between the EW and Plank mass scale. We have 
used two-loop RG equations [4T] . 


I3xs 


(dvr) 2 
1 

(47r)^ 


-1^9? + Isl - 6A - 6y} - 12A, 


Kg -|- 4:K„ 


+ 


(dvr)^ 


kI + 36A^) -I- 


(dTr)'^ 


(2.5) 


where gt refer to the SM gauge coupling, and yt denotes the top Yukawa coupling. As 
expected from Eq. (l2.5p . the correlation between the upper bounds on Hg and A* is mildly 
sensitive to s scalar mass, but strongly to each other. Our scan shows that the critical 
values are given by tig < 0.6 and A^ < 0.131. 

Of course, these bounds can be relaxed when the high energy scale above which the 
model is not effective is smaller than our reference value. 

^The vacuum is given by {H) = vl\/2 and (s) = 0 for both positive A^, and which may be violated 
when both of them are negative. However, as we will shown in Sec.3, neither negative Kg and nor negative 
As are favored by the EW vacuum stability. In this paper we choose both positive As and Kg for our 
discussion. 
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Figure 1: Upper bounds on the EW values of Kg and from the requirement of pertur- 
bativity below the Plank scale. 

3 Electroweak Vacuum Stability 

It was shown in |3T] that the EW vacuum is metastable for the observed Higgs mass 
rrih — 125 GeV, as seen from the RG equation for A. Since the sign of beta function (3x 
is always negative in the SM, A becomes negative above some critical RG scale /r*. /i* is 
approximately equal to the value h*, at which the Higgs potential is maximal. For the 
central values of top quark mass = 173.2 GeV and structure constant of QGD gauge 
coupling a = 0.1184, h* ~ 10® — 10^® GeV at the two-loop level. The situation can be 
improved to h* ~ 10^® GeV for lighter rrit within la uncertainty [12]. So, in the most of 
SM parameter space EW vacuum is only metastable. 

For EW vacuum being metastable, it suffers from large quantum fluctuation in the 
early Universe, especially during and after inflation, which leads to strong theoretic con¬ 
straint on the model parameters. 

Let us hrstly discuss the survival probability of EW vacuum during inflation. Denote 
P{h, t) the probability for Higgs held with the value h at time t, which can be determined 
via solving the Fokker-Plank equation [321 [331 [3l| . 

dP _ d 
dt dh 


Ht dP V'ih) ^ 
+ —^P 


Stt^ dh 


3Hr 


(3.1) 
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Here He refers to the Hubble parameter □, which is approximate constant during inflation. 
V' denotes derivative of Higgs potential over Higgs field. For the case h < the quantum 
fluctuation-term dominates the classical one-the last term in Eq. fl3.ip . If one ignores 
this classical term, Eq. fl3.ip can be easily solved in terms of separating variants and 
resonable boundary conditions. See [ 331133 ] for details. Then, the probability for staying 
at the domain h > other than EW vacuum at the end of inflation is given by, 

P{\h\>K) = l- r P{h,^)dh, (3.2) 

J-h, h/c 

where ~ 55 — 60 is the number of e-folds. 

logio(h*[GeV]) 



Figure 2: The value of h* at which the Higgs potential is maximal in the parameter space 
of Ks — ms- The green curves correspond to the case of SM. The black curves show the 
critical value as required by EW vacuum stability. It clearly shows that the value of h* 
can be uplifed into the region surrounded by the two black curves. 

If the probability for staying at the metastable EW vacuum is not zero as required, 
P{\ h \> h*) should be smaller than , which gives rise to, 

K > \ -—He ^ 25He ~ 2.5 X lO^^GeV. (3.3) 

V 2 TT 

^The magnitude of Hubble constant is directly related to the scalar-to-tensor ratio r, with He — 
1 X 10^® X (r/0.01)^/"' GeV. We take He = 10^® GeV as the reference point, which is a good approximation 
unless r is extremely small. 
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Obviously, this condition cannot be satisfied in the most of parameter space of SM. The 
constraint in Eq. fl3.3p can be relaxed in the context of SM by introducing a new interation 
\ H 1“^ |3ll 03] between the Higgs and Ricci scalar TZ. 

With no need of parameter ^h, this reduction can be alternatively realized in the 
HDM. It is easy to understand in terms of the new contribution to beta function 
which reads at the two-loop level, 






+ 


(dvr)^ (dvr) 


- 4Kf). 


(3.4) 


The value of h* can be uplifted for small Kg, as the sign of the new contribution in Eq. fl3.dp 
can be positive in this region. In contrast, h* may be reduced in the region of large Kg, 
where the sign of the new contribution in Eg. (13.dp is negative. The green curve in Figj2] 
clearly shows the critical value ~ 0.9. Below this critical value h* can be larger than 
SM value More interesting, as shown by the red region corresponding to the range 
Kg ~ 0.25 — 0.8, the required value = 2.5 x 10^^ GeV in Eq. fl3.3p can be satisfied. 

The effect on the vacuum stability of the singlet scalar potential due to quantum 
fluctuations is more easily understood, in comparison with the Higgs potential. Since 
(s) = 0 is the true minimal of singlet scalar potential, quantum fluctuation with magnitude 
of order 6s ~ may push s towards to an unstable vacuum from the origin. But it rapidly 
returns to the origin vacuum through classical tunning process. This understanding holds 
as long as A* is always positive between the EW and Plank scale. 

Once inflation ends, (p)reheating begins immediately. The quantum effect on the EW 
vacuum during (p) reheating can be studied in terms of the equations of motion for s 
and Higgs in principle. Unfortunately, a concrete constraint on the model parameters 
in HDM can not be derived without knowledge of reheating temperature Tre, which is 
an important physical quantity involved |3l]. To determine Tre interactions between the 
infiaton scalar and SM fields should be given explicitly. 


4 Dark Matter Phenomenology 

In this section we firstly consider the constraints arising from the DM relic density and 
direct detection limits at the LUX and XenonlOO experiments. Then we discuss the 
prospect at the XenonlT experiments. 

^For negative Kg, | | is larger in comparison with positive Kg. However, its magnitude is more 

seriously bounded above by the analysis of perturbativity. 
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4.1 Relic Density 

The Plank and WMAP 9-year data have measnred the DM relic density in high precision 

m, 

= 0.1199 ±0.0027 (4.1) 

With the assnmption that the relic density of s scalar is totally produced by the thermal 
freeze-out process, we have 

~ 0.1 pb/ {aannV) 

where (Tann is the total annihilation cross-section for s and v is the relative velocity. 



Figure 3: Contours of s scalar relic density as required by the DM projected to the 
two-parameter plane of {nis, Ks)- 

The ss annihilation mainly proceeds through Higgs-mediated process in the s channel. 
Sub-dominant process include s exchange in the t channel and annihilation into hh. So, 
(Tann IS Sensitive to the model parameters Kg and m^. 

Instead of analytic method we use the code MicrOMEGAs [38] to calculate the s scalar 
relic density with the Feynman rules generated by the package LanHEP 3.2.0 |39j . 
In FigjSjwe project the contour of relic density Vtgh? as required by Eq. fl4.ll) to the two- 
parameter plane of (m^. Kg)- As expected the required value of Kg is minimal for mg near 
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mh/^., because near this mass region the annihilation cross section is resonantly enhanced. 
Plot similar to FigJ2]is also shown, e.g., in [2T] 

Figl3] indicates that a wide mass range for s scalar is viable. However, it also shows 
that Ks bigger than 0.1 is required when rris is far away from the resonant mass region. 
For example, large Kg > 0.5 is needed for above 3 TeV, which is not favored neither 
by the EW vacuum stability nor the perturbative analysis. In the next subsection, we 
discuss the direct detection in the allowed parameter space. 



Figure 4: Spin-independent nucleon-DM scattering cross section as function of m^. Vari¬ 
ous experimental limits are also shown for comparison. 


4.2 Direct Detection 


The spin-independent nucleon-DM scattering cross section is given by. 


C^SI = 


ATTmjjnl ’ 


(4.2) 


where is the nucleon mass, fi = msm]\f/{ms + m^) is the DM-nucleon reduced mass, 
and /at ~ 0.3 [21] is the hadron matrix element. 


^ Although these plots match very well, we would like to mention that contributions to aann due to 
three- and four-body final states or QCD corrections are not included in MicrOMEGAs, with which the 
shape around W boson mass is smoother |40j . We thank the referee for reminding us this point. Since 
DM mass region near 80 GeV has been excluded by the LUX experiment, these corrections are expected 
to not affect our main results in Sec.6. 














The current best limit on usi comes from the XenonlOO [23] and LUX [25] experiments. 
In FigJUwe plot the spin-independent nucleon-DM scattering cross section as function of 
rris, with various experimental limits for comparison. The red curve corresponds to the 
DM relic density, which indicates that there are three viable regions, 

low mass region : 1 GeV < < 6 GeV, 

resonant mass region : 56 GeV < rUg < 66 GeV, 

high mass region : >185 GeV. (4.3) 

By combing the constraint from SM vacuum stability, one hnds that 

• The resonant mass region is totally excluded, which is more powerful than the 
indirect constraint from the Fermi-LAT [27] . 

• The high mass region is further reduced to more narrow mass window, which is the 
subject of Xenon IT [2l] experiment for below 4 TeV. 

• The low mass region can be consistent with the constraint from SM vacuum stability. 
However, we will show in Sec.V that the low mass region is excluded by the Higgs 
invisible decay h —)■ ss at the LHG experiment (see also some literature in Refs. [ 8 - 
18]). 

5 LHC Phenomenology 

Since the scalar singlet only couples to the Higgs directly, the signal searches for this 
scalar at the LHG mainly focus on the Higgs invisible decay h ^ ss in the low mass 
region and resonant mass region with < mh/2. Note that for scalar singlet DM there 
is no mixing effect @ between s and h, as clearly shown in Eq. fl2.ip and Eq. fl2.2p . 

Very recently, the ATLAS Gollaboration has reported the latest data about the upper 
bound on the Higgs invisible decay width [26] . 

Fi„. < 0.29 Fsm, (5.1) 

^If s scalar is not DM, one may discuss its phenomenology by introducing mixing effects. For a recent 
discussion, see [44] . 
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Figure 5: Contours of invisible decay width r(h —>■ ss) in the two-parameter plane of 
nis — Ks- The contour in black denotes the latest upper bound in Eq. fl5.ip . Region above 
this line is excluded. 

which is obviously stronger than what was used in some earlier literature. An update is 
thus meaningful. In Eq. fl5.ip Fsm — 4.15 MeV for 125 GeV Higgs mass, and the invisible 
decay width Finv is given by, 



(5.2) 


We show the decay width T[h —)■ ss) in FigJS] The contour in black denotes the 
latest upper bound in Eq. fl5.ip . and region above it is excluded. In particular, this bound 
excludes the whole low mass region and the resonant mass region with < 62.5 GeV in 
Eq. fl4.3p . In contrast, the high mass region is less constrained at the LHC in compared 
with the low mass region. 

6 Conclusions and Discussions 

In this paper we have reconsidered the simplest HDM. By combining the constraints 
arising from DM relic density, direct detection limit on spin-independent nucleon-DM 
scattering cross section from the LUX, and limit on the Higgs invisible decay from the 
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Figure 6: Mass regions consistent with various constraints. In the light of the SM vacuum 
stability, only mass region 1.1 TeV < < 2.5 TeV is viable, and it slightly reduces 

to 1.1 TeV < rus < 2.0 TeV by the additional requirement of perturbativity. For the 
explanation about h* and K^max, see previous discussions. 

LHC, we have obtained two viable mass regions (see Figj6]): 

resonant mass region : 62.5 GeV < rUs < 66 GeV, 

high mass region : rUg > 185 GeV. (6.1) 

We also show that if the constraint from the SM vacuum stability is imposed, we arrive 
at the main conclusion- the resonant mass region is totally excluded and the high mass 
region in Eq. fl6.ip is reduced to a narrow window, 

high mass region : 1.1 TeV < < 2.55 TeV, (6.2) 

as shown in FigJHl If one further imposes the constraint due to perturbativity, the mass 
range in Eq. fl6.2p is slightly reduced to, 

high mass region : 1.1 TeV < < 2.0 TeV. (6.3) 

Fortunately, as shown in FigJUthe mass region in Eq. fl6.3p can be totally detected by 
the Xenon IT experiment in the near further. Similarly to Xenon IT experiment, this 
high mass region can be also directly detected at colliders in principle. As a result of 
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small production cross section of order less than 1 fb for the s scalar at the 14-TeV LHC 
[ISj, a 5 ct discovery acquires luminosity at least of order 0(10) ab~^. 

Independently, effects on large scale structure induced by DM self-interaction is a new 
window for the detection. In compared with the abelian vector- or fermion-like DM, s 
scalar with quartic self-interaction is very distinctive. DM self-interaction has been used 
to interpret m the observations on four galaxies in the core of galaxy cluster Abell 
3827 [l6]. If the new observations indeed arise from the DM self-interaction other than 
astrophysical artifact, the high mass region is highly constrained, see, e.g., for recent 
discussion |45] . 

In variety of contexts the high mass region in Eq. fl6.2p may be modified. For exam¬ 
ple, it can be relaxed if one employs additional mechanism to stabilize the SM vacuum 
by introducing more new parameters, such as the non-minimal interaction between the 
Higgs and Ricci scalar with non-negligible coupling constant. In contrast, this high mass 
region may be alternatively strengthened if one assumes [21] that the DM relic abundance 
is partially other than totally saturated by the scalar singlet. In this case scalar s may 
play a positive role in the electroweak baryogenesis [40] . 
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